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Abstract. We have analyzed the properties of a sample 
of 33 groups and clusters of galaxies for which both opti- 
cal and X-ray data were available in the literature. This 
sample was built to examine the baryon content and to 
check for trends over a decade in temperature down to 1 
keV. 

We examine the relative contribution of galaxies and ICM 
to baryons in clusters through the gas-to-stellar mass ratio 
(Mgas/Mi,). We find that the typical stellar contribution 
to the baryonic mass is between 5 and 20%, at the virial 
radius. The ratio (Mgas/Af*) is found to be roughly inde- 
pendent of temperature. Therefore, we do not confirm the 
trend of increasing gas-to-stellar mass ratio with increas- 
ing temperature as previously claimed. 
We also determine the absolute values and the distribu- 
tion of the baryon fraction with the density contrast S with 
respect to the critical density. Virial masses are estimated 
from two different mass estimators: one based on the 
isothermal hydrostatic equation (IHE), the other based 
on scaling law models (SLM), the calibration being taken 
from numerical simulations. Comparing the two methods, 
we find that SLM lead to less dispersed baryon fractions 
over all density contrasts and that the derived mean ab- 
solute values are significantly lower than IHE mean val- 
ues: at (5 = 500, the baryon fractions (gas fractions) are 
11.5-13.4 % (10.3-12 %) and ~ 20% (17 %) respectively. 
We show that this is not due to the uncertainties on the 
outer slope f3 of the gas density profile but is rather in- 
dicating that IHE masses are less reliable. Examining the 
shape of the baryon fraction profiles we find that cluster 
baryon fractions estimated from SLM follow a scaling law. 
Moreover, we do not find any strong evidence of increasing 
baryon (gas) fraction with temperature: hotter clusters do 
not have a higher baryon fraction than colder ones, neither 
do we find the slope (3 to increase with temperature. 

The absence of clear trends between /& and Mgas/M^, 
with temperature is consistent with the similarity 
of baryon fraction profiles and suggests that non- 
gravitational processes such as galaxy feedback, necessary 
to explain the observed luminosity-temperature rela- 
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tionship, do not play a dominant role in heating the 
intra-cluster gas on the virial scale. 



1. Introduction 

Clusters of galaxies are fascinating objects because 
their observations can in principle allow one to constrain 
the parameters of the standard cosmological model. In 
particular, they are widely used as indicators of the mean 
matter density of the universe. Galaxy clusters have been 
shown to harbour very large quantities of dark matter 
since the pioneering work of Zwicky (1933), but its exact 
quantity, its spatial distribution and above all its very na- 
ture are still awaiting answers. 

Clusters are the most massive objects for which both 
the luminous baryonic mass (consisting of the X-ray emit- 
ting intracluster gas and the visible part of galaxies) and 
the total gravitating mass can be estimated. Most of- 
ten, the assumption of isothermal hydrostatic equilibrium 
(IHE) of the intra-cluster gas within the dark matter po- 
tential well is adopted to derive the total mass of clusters 
from X-ray observations, although many clusters exhibit 
obvious substructures, both in the galaxy distribution and 
in the X-ray emission morphology. 

Beyond the classical M/L ratio, clusters are at the cen- 
ter of new cosmological tests of the mean density, which 
are different in spirit and which are more global. Partly be- 
cause of this new perspective, general observational prop- 
erties of clusters have been investigated in detail in re- 
cent years. These studies were triggered by analytical ar- 
guments as well as numerical simulations which indicated 
that clusters might have similar properties in their struc- 
ture. A first means of determining the mean density from 
clusters is to use their abundance as well as their rela- 
tive evolution with redshift (Oukbir & Blanchard 1992 ; 
Bartlett 1997). A further important property of clusters 
is that their baryon fraction /b is expected to be identical 
(White et al. 1993), reflecting the universal baryonic con- 
tent of the universe. As primordial nucleosynthesis calcu- 
lations provide very strong constraints on the value of the 
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baryonic density parameter fib, determining the baryonic 
fraction in galaxy clusters allows to derive the matter den- 
sity parameter Qq = Qh/fh- This surmise, when applied 
to a set of clusters, leads to a high mean baryon fraction 
/b, of the order of 20 h^Q^^% (David et al. 1995 , here- 
after D95; White & Fabian 1995 ; Cirimele et al. 1997 ; 
Evrard 1997). Consequently, the critical value f^o = 1 is 
disfavored (as the primordial nucleosynthesis is indicative 
of fib = 0.0776 h^Q r7io/5.3 ± 7%, one obtains Qq ~ 0.4). 
White et al. (1993) have reviewed this critical issue in the 
case of the Coma cluster. 

Some caution is necessary though, since there exists an 
appreciable dispersion in the range of published baryon 
fractions. This scatter may be due to intrinsic dispersion 
in baryon fractions of different clusters. If real, it is impor- 
tant to understand the origin of such a scatter. However, 
Evrard (1997) did not find any convincing evidence for 
a significant variation in the baryon fraction from cluster 
to cluster. Such a result is in contrast with Locwenstein 
& Mushotzky (1996) and D95. These latter authors, from 
their study of ROSAT PSPC observations of a sample 
of groups and clusters of galaxies, have found a correla- 
tion between the gas fraction and the gas temperature, 
breaking the simplest self-similar picture (the different 
conclusion of Evrard could be due to the limited range 
of temperatures he used). A possible explanation for such 
variations, if real, could in principle be the development 
of a segregation between baryons and dark matter occur- 
ring during the cluster collapse, operating more efficiently 
in massive clusters. However, this mechanism has been 
shown by White et al. (1993) to be insufficient to signifi- 
cantly enhance the baryon fraction and it is therefore un- 
likely that such a phenomenon could lead to a substantial 
scatter in baryon fractions. Another possibility is that in 
poor clusters and groups, a part of the gas has been swept 
away in the shallow dark matter potential well by galac- 
tic winds, being thus less concentrated than in massive 
clusters. This scenario would also be consistent with the 
claim that the gas to stellar mass ratio increases mono- 
tonically with the temperature of the cluster (David et al. 
1990, hereafter D90). Finally, a further possibility is that 
mass estimates are not accurate and that a systematic bias 
exists with temperature. In any case, D95 derived this cor- 
relation from a very reduced set of objects (7 clusters and 
4 groups) and it would deserve further investigation based 
on a larger sample. 

As a consequence, it was one of our aims to address 
these questions with improved statistics. Moreover, in the 
baryon problem, the reliability of mass estimates is rather 
crucial and assumptions such as equilibrium and isother- 
mality may introduce systematic differences in the results 
that we wish to examine in detail. The validity of mass 
estimates has been questioned by Balland & Blanchard 
(1997). We have therefore taken the opportunity of this 
study to perform a comparison between the standard mass 
estimate based on the IHE /3-model and an alternative 



method derived from scaling arguments and numerical 
simulations including gas physics (see section 3.3), here- 
after called the scaling law model (SLM). 

In this paper, we present an analysis of a sample of 
26 galaxy clusters and 7 groups taken from the literature. 
We required that optical data were available for our ob- 
jects and searched for a precise information on the galaxy 
spatial distribution and luminosity function, on the X-ray 
temperature and on the gas density profile, in order to be 
able to build up the density and mass profiles for galaxies, 
gas and dark matter. This allows to compute properly the 
baryon fraction rather than only the gas fraction as is often 
done. This is especially important for low mass objects, in 
which the stellar component is generally believed to be 
relatively more important. Our sample comprises clusters 
with temperature from 1 to 14 keV, and therefore allows 
us to investigate several interesting quantities beyond gas 
and baryon fraction, like the mass to light ratio and the 
ratio of galaxy baryonic mass to gas mass (possibly pro- 
viding important constraints on galaxy formation), over a 
wide range of temperatures. All the data used here come 
from the literature, with the exception of Abell 665, for 
which we have analysed an archival ROSAT image to ob- 
tain the gas density profile. In fact, this cluster has already 
been studied from Einstein data by two teams (Durret et 
al. 1994 ; Hughes & Tanaka 1992), finding in each case 
a surprisingly very high gas fraction (respectively ~ 50% 
and 33%). We will see this cluster provides a striking ex- 
ample of the differences in mass determinations. 

The sample is presented in section 2. The methods to 
compute the various quantities for each cluster in the sam- 
ple is presented in section 3 and the results are presented 
in Section 4. In section 5, we examine the trend with tem- 
perature for several quantities. 

In all the present study, we assumed a Hq = 
50 kms~^ Mpc~^ and qo = 0.5 cosmology. 

2. The sample 

We looked in the literature for objects studied thor- 
oughly enough to allow us to compute the baryonic mass 
in galaxies, the mass in gas and in dark matter at any 
radius. This is quite not a refinement, since both the 
baryon fraction and the galaxy to gas mass ratio can 
vary very rapidly with radius, as will be seen in the next 
section. We therefore needed detailed information, which 
drastically reduced the possible number of objects that 
could be included in the sample. When a same object 
was studied by several teams, wc applied straightforward 
selection criteria : for spatial X-ray data, for instance, 
we systematically prefer ROSAT observations, because of 
its improved spatial resolution and sensitivity, whereas 
for X-ray temperatures, Ginga and ASCA satellites were 
preferred to Einstein MFC, most of which come from 
the catalogue of David et al. (1993). Recently, it has 
been noticed that cluster luminosities and temperatures 
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might change noticeably when the central cooling flow 
emission is removed (Markevitch 1998 ; Arnaud & Evrard 
1999). It is not clear which temperatures are to be used 
(especially when using a mass-tcmpcraturc relationship 
derived from numerical simulations). In order to keep 
our sample as homogenous as possible, we did not use 
cooling flow-corrcctcd temperatures which arc not always 
available. Furthermore Markevitch (1998) found that 
temperatures corrected for central emission are in the 
mean 3% larger, which will be of weak consequence in 
our average quantities. However, our treatment of the 
uncertainties on temperatures leads to large error bars 
when a large dispersion in measured temperatures exists 
(see Table 1), as for instance in the presence of strong 
cooling flow. 

In some cases, optical data may be very uncertain 
because of projection effects and magnitude limitations, 
especially for groups whose galaxy membership is some- 
times tricky to establish. However, we tried to identify 
objects for which data are reasonably reliable and we 
derived mean dynamical quantities for this sub-sample 
as well. Finally, it must be emphasised that the X-ray 
limiting radius at which baryon fractions are estimated is 
a crucial parameter, since both the galactic mass derived 
from a King profile and the X-ray gas mass given by the 
Hubble-King model diverge respectively for e < 1 and 
/3 < 1 (the definition is given in 3), requiring that they be 
truncated. It is also important that the baryon fractions 
of different clusters be computed at an equivalent scale 
in order to test the scaling hypothesis and if statistical 
conclusions are to be brought out from them, i.e. that 
wc use the radius containing the same overdensity, while 
information is actually available only up to the X-ray 
limiting radius i?xiim which primarily depends on the 
characteristics of the observations (detector sensitivity, 
integration time...). 

X-ray and optical data are summarised in Tables 
1 and 2, using a Hubble constant /150 — 1. Notes on 
clusters which required a special treatment due to an 
incompleteness of data can be found at the end. Optical 
luminosities are given in the blue band. When the blue 
luminosity was not available, wc used the following colors, 
corresponding to standard values for elliptical galaxies : 
B-V = 0.97, V-F = 0.76, r-F = 0.58 (Schneider et al. 
1983) and R = F (Lugger 1989). 



2.1. The case of Abell 665 

This cluster is one for which large baryon fraction esti- 
mates have been published in the literature. As these are 
surprisingly high, we have found interesting to re-analyse 
this cluster using a ROSAT archival image and the cali- 
bration routines of Snowden et al. (1994). We found that 
the gas surface brightness profile is well fitted by a Hubble- 
King law, and the X-ray emission can be traced out to a 




radius (arcmin) 

Fig. 1. Surface brightness profile of the intracluster gas 
in A665, in the ROSAT bands R4 to R7 (0.44 to 2 keV). 
Points represented by an empty circle have been excluded 
from the fit because at these radii some background or 
foreground X-ray sources appear in the map. The dashed 
line is the fitted background level. 



very large radius. The background level, which has been 

fitted together with the other parameters, is estimated 
with comfortable confidence. Spherical symetry was as- 
sumed to derive the surface brightness profile in 0.5 arcmin 
wide annuli, although the X-ray map shows significant 
departure from sphericity ; however, the effect of ellip- 
ticity on derived masses is known to be negligible (Buote 
& Canizares 1996). The central electron volume density 
rieo was computed by matching the theoretical count rate 
with the 0.547 counts s~^ collected within a 30 arcmin ra- 
dius (after subtraction of the background), which amounts 
to solving: 



ntr. aE D T 



g{T,E) e-^ g-'^^^)^" A{E) 



E 



dE 



2\(-3/3+i) 



Ode, 



(1) 



with aE = 1.02 lO"'^^ SI and the angular distance 
D ~ 812 Mpc, and where A[E) stands for the en- 
ergy dependence of the transmission efficiency. The re- 
sults of this analysis are the following (for the bands R4 
to R7 of ROSAT) : /3 = 0.763 ± 0.023, 6^ = (112 ± 
5) arcsec (which corresponds to 0.44 Mpc at the distance 
of A665), neo = (2.85 ± 0.38) 10"^ electrons cm'^ and 
-Rxiim = 10 arcmin (= 2.36 Mpc) with a central surface 
brightness Sq = (3.53 ± 0.26) 10~^ counts s~^ arcmin"^ 
and a background surface brightness 6 = (2.2 ± 
0.6) 10~^ counts s~^ arcmin^^. We used the gas tempera- 
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ture and foreground absorbing hydrogen column density 
given by Hughes & Tanaka (1992) from their Ginga analy- 
sis, together with the formula of Mewe et al. (1986) for the 
Gaunt factor and that of Morrison & McCammon (1983) 
for the interstellar absorption cross section. At -Rxum, 
the inferred gas mass is Mgas = (2.46 ± 0.76) IO^Mq 
(at, which is similar to the values found by Durrct ct 
al. (1994) and Hughes & Tanaka (1992). The hydrostatic 
mass is Mhydro = (1-60 ± 0.24) lO^^M©. Our mass es- 
timate from NFW's dark matter profile, computed with 
the EMN normalisation (sec section below) is Mjark = 
(1.41 ± 0.25) IO^^Mq, and the resulting total mass is 
MsLM = (1-68 ± 0.33) IQI^Mq. The baryon fraction 
amounts to respectively (16.3 ± 7.5)% and (15.6 ± 6.4)%. 
Hence A665 is a quite ordinary rich cluster whose baryon 
fraction seems reasonable if compared to previous values. 
We have also compared our gas mass estimates for the 
whole sample with other published analyses and found 
good agreement while the main differences are on /{,, com- 
ing from the estimation of total masses as will be discussed 
in 5. 

3. Analysis methods 

3.1. The stellar mass profile 

The stellar matter content can be computed at any 
radius from the cluster center using the projected number 
density profile of galaxies, their luminosity function and a 
mass to light ratio for the stellar population calibrated on 
the observation of nearby galaxies. Most often, the density 
profile is fitted by the common King form : 



CTgal(r) = O-o (l + (^' 



with 



(2) 



where rc is the galactic core radius. The case e = 1 is an 
approximation to the isothermal sphere, in which galaxies 
have reached their equilibrium distribution. The advan- 
tage of such a model is that the volume density is obtained 
by an analytical deprojection. However, de Vaucouleurs 
profiles, which are much steeper in the cluster core, pro- 
vide a better approximation to the real distribution (Rhee 
& Latour 1991 ; Cirimele et al. 1997), at the same time 
leading to a finite total number of galaxies : 

CTgai(r) = (JO exp (- {■^y'^ . (3) 

This sort of profile was deprojected using the formula : 



v{r) 



1 d_ 
2'Kr dr 



r+oo 



dp^ 



(4) 



RXJ 1340.6+4018 consisting of only one giant elliptical 
galaxy, for which we used Mt^/L^ = 8.5 /isqMq/Lb©) 
is M*/iB = 3.2 /isqMq/Lbo, obtained by White et al. 
(1993) by averaging over the Coma luminosity function 
the M*/I/ ratio from van der Marel (1991) given as a 
function of luminosity for bright ellipticals. Then, using 
the Schechter luminosity function : 



n{L)dL = N* (^)-%-r-d(^). 



(5) 



the luminosity emitted by a shell of thickness dr and sit- 
uated at the radius r writes as : 
dL{r) = Ltot X (47rr2 i'gai(r)dr) / A^(> Lum) 
where N{> L) = N* V{1 - a,L/L*) is the total number 
of galaxies brighter than L, Liim being the limiting lumi- 
nosity of the observations, and Ltot = N*L*T{2 — a). The 
stellar mass enclosed in a sphere of radius R can eventually 
be written as : 



M,(fl) = — / 47rr2 z/gai(r)rfr. 

L r(l - a, -^) Jo 



(6) 



When no parameters for the luminosity function were 
found in the literature, we adopted the standard ones 



(Schechter 1975) : a = 1.25 and 



-21.9 -h 5 log /i5o. 



As a few clusters observed in X-rays do not have any 
available spatial galaxy distribution (or with too poor 
statistics), but only either a luminosity profile or even 
several total luminosities given at different radii, we then 
assumed a King profile and fitted the few points by the 
resulting integrated luminosity profile : 



L{< H) = Lo 



by varying simultaneously Lo and Tc . In addition to those 
cases, RXJ 1340.6-1-4018 was treated in a special way : we 
deprojected a de Vaucouleurs luminosity profile (Ponman 
et al. 1994). 

3.2. The X-ray gas mass profi,le 

In their pioneering work, Cavaliere & Fusco-Femiano 
(1976) have shown under the isothermality assumption, 
that the X-ray gas profile is described by: 



(8) 



Pgas(r-) = Po (l+ (t^i 



which translates to the observed X-ray surface bright- 
ness with the following simple analytical form (the so- 
called /3- model): 



p being the projected distance to the cluster centre and S{9) = So\\-\- (j-) 
r the true distance. Because this deprojection is numer- 
ically unstable, we computed it by assuming u{p) to be 
constant inside a grid step and then integrating analyt- 
ically the denominator. The mass to light ratio applied 
to all clusters and groups (but the supposed fossil group 



(9) 



the slope j3 and the core radius Tcx, which are interdepen- 
dant in their adjustment to the surface brightness, being 

generally found to range between 0.5 and 0.8 and between 
100 and 400 kpc respectively. Very often, central regions 
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Fig. 2. Mass and mass ratio profiles for a few objects. The meaning of line styles is as follows : Left panels : thick 
line : SLM mass ; thin line : IHE mass ; dashes : gas mass ; dot-dashed line : stellar mass. Right panels : thick lines : 

baryon fraction (continuous), gas fraction (dashed) and stellar to total mass ratio (dot-dashed) in the SLM case (with 
EMN calibration) ; thin lines : same quantities for the IHE model ; three-dots-dash : stellar to gas mass ratio. 
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Fig. 3. Profiles of the baryon fraction and gas fraction as a function of mean overdensity for objects with the most 

reliable data. Left panels show these profiles in the case of the hydrostatic assumption and right panels for mass 
estimates derived from NFW's dark matter profile, with EMN normalisation. Groups (Tx < 2 keV) are represented 
with dotted lines, cool clusters (Tx < 5 keV) with dashed lines and hot clusters with continuous lines. The group with 
a very steeply rising baryon fraction in the IHE case has /? = 0.31. 



of clusters have to be excluded from the fit, due to cool- 
ing flows resulting in an emission excess. The gas mass 
can be inferred accurately from the knowledge of So, /3 
and 9c- Uncertainties in the gas mass are small in gen- 
eral, as long as it is computed inside a radius at which the 
emission is detected. The relationship between the elec- 
tron number density and the gas mass density used here 
is po = 1-136 mp rieo (assuming a helium mass fraction of 
24% and neglecting metals). 



3.3. The binding mass profile 

Mass estimation is certainly the most critical aspect 
of recent studies of the baryonic fraction in clusters. 
Clarifying this issue is one important aspect of this paper. 
We derived the gravitational mass in two ways : 



The hydrostatic isothermal j3-model : 



First, we used the standard IHE assumption which, 
using spherical symmetry, translates into the mass profile : 



Mtot(r) 



Gunip 
3k 



Txr 



dln/)gas(r) 
cflnr 



/3Txr(l+(^)-^)"\ 



(10) 



The total mass thus depends linearly on both (3 and Tx. 
Hence, if the slope of the gas density is poorly determined 
(and this is the case if the instrumental sensitivity is too 
low to achieve a good signal to noise ratio in the outer 
parts of the cluster), it will have a drastic influence on 
the derived mass. This mass profile results in the density 
profile : 



Ptot(r') oc \ 



(11) 



and in a flat density at the cluster centre. 

The isothermality assumption can be questioned since 

Markevitch et al. (1998) foimd evidence for strong 
temperature gradients in clusters, which may lead to 
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IHE mass estimates smaller by 30% (Markcvitch 1998). 
Although the reality of these gradients has recently been 
questioned (Irwin et al., 1999; White 2000). 

• The universal density profile : 

An alternative approach is to use the universal dark mat- 
ter density profile of Navarro et al. (1995, hereafter NFW) 
derived from their numerical simulations : 



0.25 



Pdark(r) 



1500 



r (5r + r2oo)^ 



(12) 



where r2oo stands for the radius from the cluster center 
where the mean enclosed overdensity equals 200 (this is 
the virial radius) and is the critical density. It varies as 
near the centre, being thus much steeper than in the 
hydrostatic case ; NFW claim this behaviour fits their high 
resolution simulations better than a flat profile. Further- 
more, contrary to the /?-model , the dark matter density 
profile obtained by NFW is independent on the shape of 
the gas density distribution. This will introduce a further 
difference. The normalisation of the scaling laws ensures a 
relationship between temperature, virial radius and virial 
mass. Here, this normalisation is taken from numerical 
simulations. Different values have been published in the 
literature (see for instance Evrard, 1997 ; Evrard et al., 
1996, EMN hereafter ; Pen, 1998 ; Bryan and Norman, 
1998, BN hereafter). Frenk et al. (1999) investigated the 
formation of the same cluster with various hydrodynami- 
cal numerical simulations. They found a small dispersion 
in the mass-temperature relation : the rms scatter a is 
found to be of « 5%, EMN and BN lying at the edges 
of the values found, representing a 4 tj difference. EMN 
provide a scaling law between r5oo and Tx : 

r-500 = 2.48 (^) ^ /iso^ (1 + z)-i Mpc (13) 

(in terms of comoving radius) which was used here to com- 
pute r2oo, writing : 



^(r-50o) = 500 



(l^^ioo)"' / 

Jo 



Airr 



2 P(0 

Pc 



dr 



180/ X"^ [ ln(l + 5X) - 



5X 



T+5XJ (14) 

with X = r5oo/''2oo5 where / = 1739/1500 is a corrective 
factor to transform dark matter mass into total mass, so 
that 5(r2oo) is really equal to 200. Solving this equation 
gives X = 0.66. The relation at 2; = between virial mass 
and temperature can then be written as : 



Tx = 4.73Mi5(r2oo)5 keV. 



(15) 



BN did provide the following constant of normalisation : 
Tx = 3.84Mi5(r2oo)^keV. (16) 

This difference is quite significant : it does correspond to 

a virial mass 40% higher. Using this normalisation will 
obviously significantly change the inferred gas fraction. 




Fig. 4. Average profiles for all clusters (groups included) 
with the most reliable data. Top : baryon fraction (contin- 
uous lines) and gas fraction (dashes) in the case of SLM 
mass estimates with the EMN normalisation (thick lines) 
and hydrostatic masses (thin lines). Middle : mass to lu- 
minosity ratio for the whole sample (continuous lines), for 
King galaxy profiles only (dots) and for de Vaucouleurs 
profiles only (dash-dots), with the same convention as pre- 
viously. Bottom : stellar mass to gas mass ratio, with the 
same line styles as for M/L-q- 
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4. Results 

4-.1. Distribution of the various components 

The results from our analysis for each cluster in our 
sample are summarised in Tables 3 to 6 : in Table 3 and 
Table 5 mass estimates at r2oo are derived from the SLM 
with the Tx - My calibrations respectively given by EMN 
and BN (used to compute r2oo)- In Table 4 and Table 6 are 
summarised mean dynamical quantities over the sample 
at three different overdensities, r2oo, ''soo and r2ooo- The 
same quantities are also given with mass estimates from 
the IHE model. 

4.1.1. The binding mass 

Mass profiles of the various components for a few clus- 
ters are displayed in figure 0, together with mass ratio 
profiles (right side). Figure gshows baryon and gas frac- 
tion profiles for the whole sample. Quantities are plotted 
against the mean enclosed contrast density, which is the 
natural variable in the scaling model. A clear feature aris- 
ing from Figure ^ concerns the different behaviours of hy- 
drostatic masses and total masses deduced from NFW's 
dark matter profile, normalised by the EMN Tx - My re- 
lationship : NFW profiles are more centrally concentrated, 
as could be foreseen from Equ.ll and 12, a property which 
is in agreement with the density profile of clusters inferred 
from lensing (Hammer 1991, Tyson et al. 1990). In the 
outer part, when the contrast density is smaller than a 
few 10^, the shapes of the density profiles are quite similar, 
although some difference in the amplitude exists. In fact, 
profiles calibrated from the EMN Tx-My relation tend to 
be systematically more massive than with the isothermal 
hydrostatic model, with a significant dispersion. The last 
column of Tables 3 and 5 gives the ratio between masses 
computed with both methods. The mean of masses esti- 
mated by the IHE /3-model is significantly smaller than 
SLM masses (at i?xiim): Mihe/Mslm = 0.80 ±0.03 with 
EMN normalisation, and Mihe/Mslm = 0.67 ± 0.026 
with BN normalisation. Clearly, such a difference will 
translate into the baryon fraction estimates. 

4.1.2. The X-ray gas 

Second, the distribution of gas is more spread out than 
that of dark matter, which results in steadily rising baryon 
fractions with radius (Fig. as was already pointed out 
by numerous teams, among which Durret et al. (1994) 
and D95. NFW also recover this trend in their simula- 
tions. This fact makes the choice of the limiting radius an 
important matter. In particular, extrapolating masses to 
the virial radius (which is reached by the gas emission in 
only five clusters among our sample) could be very unsafe, 
especially for cool clusters, the gas of the most extended 
of our objects with Tx < 5 keV being detected only out to 
S = 500. 



4.1.3. Mass to light ratio 

The derived mean mass to blue luminosity ratio is 
shown in figure I As it can be seen, M/Lb remains re- 
markably constant from S ~ 5000 to the outer parts of 
clusters, in the case of total masses derived from SLM 
as well as that of hydrostatic masses. Thus, the widely 
spread assumption that light traces mass is confirmed, at 
least at r > r^- The influence of the choice of de Vau- 
couleurs galaxy density proflles as compared to King pro- 
files is also clearly highlighted. In fact, in the core, dark 
matter is normally much more concentrated than galaxies, 
but using a dc Vaucoulcurs distribution, it turns out that 
the concentration factor is considerably lowered and even 
reversed in the case of hydrostatic masses. Mixing the two 
shapes of galaxy distribution in our sample, the result is 
an intermediate behaviour. 

4-. 2. The baryon fraction 

We find that inside a same object, the gas and baryon 
fractions increase from the center to outer shells (Fig. || 
and Fig. p, refiecting the fact that the distribution of 
gas is fiatter than that of dark matter, a trend similar to 
what is found by D95. Secondly, an interesting feature can 
be noted from figure |^: the baryon fraction profiles versus 
density contrast are remarkably similar and seem to follow 
a regular behaviour, consistent with an universal baryon 
fraction shape, even in the central part (although with 
a larger dispersion). This behaviour appears more clearly 
when one is using the SLM model. This result is consis- 
tent with the baryon fraction following a scaling law as it 
has been already found for the emissivity profiles (Neu- 
mann & Arnaud, 1999) and gas profiles (Vikhlinin et al., 
1999). Thirdly, the comparison of the graphs of figure || 
shows that the baryon fraction fsLM estimated from the 
NFW profile normalised with the EMN Tx ~ My relation- 
ship is less dispersed at all contrast densities. This effect is 
asymmetric : the high baryon fractions fiHE found with 
the IHE method disappear. The fact that fsLM appears 
less dispersed has already been found by Evrard (1997) 
and Arnaud & Evrard (1999). However, our work indi- 
cates that this feature exists at any radius. We also plot 
in Figure |^ the histograms of baryon fractions as derived 
from both the IHE and SLM methods at the virial radius 
r2oo but also at r2ooo chosen because each object of the 
sample is detected in X-rays at least out to i5 ~ 2000. 
The comparison of the two indeed provides an evidence 
for SLM masses to lead to more tightened baryon frac- 
tions than hydrostatic masses. At the virial radius, we 
found that the intrinsic dispersion is 50% with the IHE 
and 20% with SLM. This bears an important consequence 
for the interpretation of mass estimates as well as the in- 
terpretation of the baryon fraction. Clearly the fact that 
the baryon fraction is less dispersed in the SLM even in 
the central regions shows that this mass estimate is safer 
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and that the IHE method provides less accurate mass es- 
timates even in the central region where hydrostatic equi- 
librium is expected to hold. 

4.2.1. Stellar to gas mass ratio Af*/Afgas 

Also shown in figure ||is the mean Af»/Mgas ratio as a 
function of overdensity, slowly going down after the central 
part. The galaxy density is indeed steeper than that of gas, 
decreasing in with e = 1 instead of for a typical 
value of /3 = 0.66, and the situation is even worse when a 
de Vaucouleurs profile is used for the galaxy distribution. 
Again, the latter contributes in a large amount to the steep 
decrease in the central regions, whereas there the ratio is 
flat with King profiles. 

4-.3. Numerical results 

Average numerical results are presented in Tables 4 
and 6. It is found that the mean baryon fraction using the 
SLM with the Tx - My normahzation of EMN is 13.4% 
and the gas fraction 11.5% at r^oo to be compared with 
hydrostatic results : respectively 19.2 and 17.0%. As ex- 
pected, the two methods of mass estimation lead to differ- 
ent baryon (gas) fractions. This difference is not negligible 
(« 40%) and is mainly due as already noted, to the differ- 
ence between the IHE mass and the SLM mass. The IHE 
mass can be 50 to 60% lower with respect to the SLM 
mass (this is the case, for instance, of the groups HCG 62, 
NGC 2300 and NGC 4261).This difference between f^^'^' 
and fl^^: increases when using the Ty^-My normaliza- 
tion of BN (the mean baryon fraction being then 11.5 and 
the mean gas fraction 10.3%). Cirimele et al. (1997) found 
/b = 23% for their 12 clusters included in our sample (and 
20% excluding A76), instead of our result of 19% (and 



16%) using their parameters and the same hydrostatic /3- 
model and their limiting radius (they choose a uniform 
^xiim = 1-5 Mpc ) and of 13% using the SLM method. 
The disagreement is due to the adopted stellar mass to 
light ratio {M / L-q = 10.7 instead of our 3.2 value). From 
the results of D95, it comes out that their 7 clusters also 
have a mean baryon fraction of fb ~ 23%. Thus, this is 
a confirmation of the divergence between hydrostatic (3- 
model mass estimates and SLM's masses. From a sample 
of 26 clusters among which 7 hot and 3 cool clusters are in 
our sample, Arnaud & Evrard (1999) have made a similar 
analysis and derived in the frame of simulation-calibrated 
virial masses a mean gas fraction at 5 = 500 of ~ 14% in 
rhough agreement with our value of 12%. If the compari- 
son is restricted to hot cluster subsamples, the agreement 
is as good (they found 16 % to be compared with our 
14%), and also at r2oo- A somewhat higher gas fraction 
{fg « 17%) has been obtained recently by Mohr et al. 
(1999), as compared to ours, which is probably due to the 
difference in the normalization of the Tx — My relation- 
ship. 

Another output from the present study is the mean to- 
tal mass to blue luminosity ratio M/Lb — 270 at r2ooo 
(the hydrostatic assumption leading to M/L-q — 200), 
groups and clusters of all temperatures put together. How- 
ever, when looking in more detail at the three classes of 
groups (with Tx < 2 keV) , cool clusters (Tx < 5 keV) and 
hot clusters, M / L-g, (at r2ooo) goes from 200 to 270 and 
340 respectively, with similar statistics (7 groups, 10 cool 
clusters and 8 hot clusters). Hence, we disagree with D95 
who claim that the mass to light ratio is roughly constant 
from groups to rich clusters (using the group NGC 5044 
which also belongs to our sample with M / L-g, — 160, 2 
cool clusters and 4 hot clusters, 3 of which are also in 
common with ours). It is worth noticing that 2 of the 3 
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clusters in common have a low M/Lb in our analysis : 150 
for A85 and 170 for A2063. This conclusion holds what- 
ever the limiting radius : there is a factor of 1.7, 1.9 and 
1.9 respectively between groups and hot clusters when ex- 
amining M/Lb out to r2ooo, rsoo or r2oo- 

As to the mean gas to stellar mass ratio A/gas /A/,, its 
values are summarized in Tables 4 and 6. We have com- 
puted this quantity to estimate the stellar contribution 
to the baryon fraction and to investigate any correlation 
with temperature, which will be discussed in the next sec- 
tion. Let us simply mention that our value for groups at 
^"2000 is in good agreement with the mean value ~ 5 of 
Deir Antonio et al. (1995) for 4 poor clusters, after cor- 
recting for the different M*/L they have used. 

5. Correlations of the baryon population 
properties with temperature 

In order to properly understand the baryon fraction in 
clusters it is necessary to understand what are the rela- 
tive contribution of the gas and stellar components. Sev- 
eral previous studies found that the stellar component is 
more dominant in low temperature systems, the lower gas 
content of small clusters being possibly due to feedback 
processes. Our sample, being large and covering tempera- 
tures from 1 to 14 keV, allows us to study these questions 
in detail. 

5.1. The Mgas/Af* -7x correlation 

In this section, we examine a possible correlation be- 
tween the X-ray gas temperature and the ratio of gas mass 
to stellar mass, Mgas/M*, at various radii. A strong corre- 
lation has been previously found by D90: from the analysis 
of twelve groups and clusters with temperatures ranging 
from 1 to 9 keV, D90 found that this ratio varies by more 
than a factor of five from groups to rich clusters. An in- 
crease of Mgas/A/* with cluster richness has also been re- 
ported by Arnaud et al. (1992). 

This trend has been interpreted as due to the galaxy for- 
mation being less efficient in hot clusters than in colder 
systems. D90 suggest that the scenario for structure for- 
mation of hierarchical clustering, in which large structures 
form after little ones by successive mergers, is adequate to 
explain their result : in fact, as mergers go on, the intra- 
cluster gas is progressively heated by shocks to higher and 
higher temperatures, as the size of the structures involved 
increases ; the higher Tx, the more difficult it becomes for 
the gas to collapse and to form new galaxies. Hence, after 
some time, further galaxy formation would be prevented 
in hot clusters, producing an anti-bias. 

In figure^ we have plotted Mgas/M* against temper- 
ature at the radii r2ooo and r2oo- As can been seen, mean 
figures for groups, cool clusters and hot clusters seem to 
show the sequence observed by D90, but in a less pro- 
nounced way: we find that cool clusters have Mgas/M* 



which is ~ 3 times smaller than for hotter ones instead 
of a factor of > 5 in D90. Moreover, this apparent se- 
quence weakens when plotted at ^200: Mgg^/M^ is only 
twice smaller for groups than for hot clusters. 

It should be kept in mind that in figure ^, we adopted 
a constant galactic mass to luminosity ratio for clusters 
and groups, whereas it is expected to be lower for late 
type galaxies than for E-SO. As morphological segregation 
tends to raise the fraction of early type galaxies in rich 
clusters, taking into account this variation of M^/L with 
morphological type would in fact flatten further the ob- 
served correlation between Mgas/M* and Tx, as would do 
taking into account the difference in galactic output from 
groups to clusters. We conclude that our sample does not 
show a strong evidence, if any, of increasing Mgas/M* with 
Tx as previously found by D90. 

5.2. The /b - Tx correlation 

Our sample spans an unprecedented wide range of tem- 
perature, allowing to test the somewhat puzzling evidence 
that cool clusters have lower mean gas fraction than hot 
clusters. This trend has been first reported by D95 and 
seems to be confirmed (Arnaud & Evrard 1999). A modest 
increase of the gas fraction with Tx has also been reported 
by Mohr et al. (1999). Such a trend is unexpected in a self- 
similar cluster evolution, /gas and /b at a given overdensity 
being expected to be constant, but would be naturally ex- 
plained by non gravitational processes such as galaxy feed- 
back (for instance, early supernovae-driven galactic out- 
flows) , able to heat the intergalactic gas enough to make it 
expand out (Metzler & Evrard 1994, 1997; Ponman et al. 
1999). This is achieved more easily in shallower potential 
wells like those of groups, which could even experience 
substantial gas expulsion, thus reducing their gas frac- 
tions. Such scenarii are necessary to explain the Lx — Tx 
relation (Cavaliere et al., 1997). 

In order to examine this issue, we plot in figure ^ the 
baryon fraction versus the temperature at different radii: 
Rxiim, ^"2000 and r2oo. Error bars were estimated by con- 
sidering uncertainties on the temperature and on metallic- 
ity for groups. Uncertainties on X-ray emission are small 
and leads to tiny uncertainties on the gas mass in the ob- 
served range (i? < Rx), while in the outer part, where ob- 
servations are lacking, robust estimates of the uncertain- 
ties cannot be obtained, given that these uncertainties are 
systematic in nature. In the case of groups, metallicity un- 
certainty can lead to significant uncertainties on gas mass, 
and was therefore taken into account. As it can be seen, we 
do observe no obvious trend with Tx- The data are more 
consistent with /b being constant and this whatever the 
mass estimator used. Although a weak tendency could be 
seen (in the frame of SLM masses), it appears swamped 
in the high dispersion affecting objects of a same temper- 
ature. Therefore we do not confirm the trend of increasing 
fb with Tx (or size) as previously found by D95. This is a 
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rather robust conclusion as our sample covers a wide range 
of temperature, from 1 to 14 keV. This result is consistent 
with the similarity of the baryon fraction profiles wc found 
(Fig. ^) and the absence of trend of Mgas/A^* with Tx in- 
dicating that non-gravitational processes such as galactic 
feedback are not dominant in determining the large scale 
structure of the intracluster medium. 



5. 3. The fh~P correlation 

Analysing the baryon fraction versus temperature may 
hide or refiect some correlation which are present among 
other parameters. Of special interest is to check whether 
a correlation with /3 exists. 

We first searched for a trend between /3 and the tem- 
perature. Previous studies have shown that low temper- 
ature systems exhibit a more extended ICM distribution 
(low f3 values) than hotter ones (Arnaud & Evrard 1999). 
From figure || we can see that no clear trend of increas- 
ing (3 with Tx is found. Although smaller (3 are found at 
the cool side, this might be due to a larger dispersion in 
P for the smallest potentials. We note that our result is 
consistent with the recent analysis of Mohr et al. (1999). 

We have also examined the way the baryon fraction 
varies with f3 (Fig. H). The baryon fraction derived from 
the hydrostatic (3 model, fl^^ , does not vary with (3 in 
an obvious way: if anything it decreases with increasing [3 
(while no obvious correlation is found with Tx, see Fig. ^. 
Such a trend if real would be unexpected. Using the SLM 
mass estimates, the baryon fraction is much more constant 
and less dispersed, even at a fixed (3 (for /3 ~ 0.5 — 0.7 the 
dispersion on fb/ fb is 0.23 with SLM estimates while it 
is 0.31 with the IHE). The fact that is constant with 
(3 again differs from what one would expect if reheating 



would have a dominant role in redistributing the gas inside 
clusters. 



5.4. Implications on mass estimates 

As we have seen (section 4.2) the baryon fraction esti- 
mated with the SLM method is less dispersed than with 
the IHE method. This effect has been noticed previously 
(Evrard, 1997) and has been interpreted as due to the ob- 
servational uncertainties in the estimation of (3. The mass 
estimates (at some radius R) can be written as: 



Ml = aiTR 

and for the IHE model 

M2 - a2l3TR. 



(17) 



(18) 



The fact that baryon fractions estimated with M2 are 
more dispersed can be understood just because of the ex- 
tra dispersion introduced by f3 (EMN ; Arnaud & Evrard, 
1999). For this to be due to the sole errors in the mea- 
surement of P, it would imply that the dispersion in the 
measurements dominates the intrinsic dispersion, result- 
ing in a tight correlation between (3 and /b, which is not 
obvious from Fig. |^: most clusters have a /3 in the narrow 
range 0.5 — 0.7, and the sample restricted to this range 
shows a larger dispersion for the baryon fraction computed 
with the IHE. Therefore, we conclude that the large dis- 
persion observed in the baryon fractions estimated from 
the hydrostatic /3 model is intrinsic to the method itself 
leading to less reliable mass estimates, rather than to the 
uncertainty on /3 measurements. 

6. Discussion and conclusions 

We have analysed a sample of 33 galaxy clusters and 
groups covering a wide range of temperatures. For all 
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clusters, X-ray and optical data were gathered from the ponents of X-ray clusters. Mass estimates were derived 

literature (except for Abell 665 whose ROSAT PSPC from two different methods : first we have followed the 

data have been reanalysed by us). This has allowed us standard hydrostatic isothermal equation (IHE method), 

to investigate the structure of the various baryonic com- secondly we have estimated the virial mass and mass 
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profile by using the universal dark matter profile of NFW 
in which the virial radius is deduced from the sealing 
relation argument (SLM method), the normalisation 
constant being taken from EMN and from BN. 

We find that virial masses (i.e. masses enclosed inside 
a fixed contrast density radius) are systematically and 
significantly lower when one is using the hydrostatic 
isothermal equation. After this paper was submitted, we 
have been aware of a recent similar study by Nevalainen et 
al. (2000) who found that taking into account temperature 
profiles exacerbatc^s this difference, as inferred masses arc 
then smaller. Examination of the baryon fraction versus 
contrast density has shown that the baryon fraction is 
more dispersed using th(^ IHE. We have shown that this is 
not due to uncertainties on the /3 measurement but rather 
reflects the fact that the IHE method does not provide as 
reliable a mass estimate as the SLM, neither in the inner 
parts nor in the outer regions. Moreover the tightening 
of fh^^ profiles supports the idea that baryon profiles 
in clusters do have a rather regular structure, i.e. that 
gas distribution is nearly self-similar which is consistent 
with the recent studies by Vikhlinin et al.(1999) and by 
Neumann & Arnaud (1999) who found an evidence of 
regularity in gas density profiles. However, when plugging 
their mean standard density profile into the hydrostatic 
equation, these last authors found a mean total mass 
profile which is different from the NFW profile (their 
mass profile is lower than the one derived from numerical 
simulations). They have used the hydrostatic isothermal 
equation to estimate their total mass which probably 
explains this discrepancy. 

Our mean gas fraction at the virial radius r2oo, using 
the SLM, is found to be in the range 12.6-14.6 % (for 
h = 0.5), in rough agreement with Arnaud and Evrard 



(1999) when EMN normalization is used. The mean 
baryonic fraction is Z^^^*^ « 13.7 — 16.4. It is important 
to emphasize that our analysis shows that a larger baryon 
fraction could be obtained when the sole hydrostatic 
equation is used, but it is reasonable to think that this 
is an overestimation due to the mass estimator itself. 
Our analysis is consistent with an intrinsic dispersion of 
20% in baryon fractions (but this could be due to some 
systematics) , which means that our mean baryon fraction 
is uncertain by less than 0.01. 

As the observed luminosity-temperature does not fol- 
low the simple scaling expected from self-similarity, it 
is likely that non-gravitational heating such as galactic 
winds, additional energy input by Type 11 supernovae, 
play an important role in the physics of the X-ray gas 
and may result in inflating the gas distribution. Metzler 
& Evrard (1997) have studied this possibility and found 
that this is achieved more easily in low temperature clus- 
ters (shallow potential wells). The consequence of such 
an effect is an increasing gas fraction outwards within a 
cluster and a decreasing gas fraction with decreasing tem- 
perature. This effect is expected to be more pronounced 
in groups and cold clusters. From our sample, we do not 
observe such a trend of the baryon (gas) fraction with Tx 
whatever the method we use, suggesting that non gravita- 
tional heating is not playing a dominant role on the scale 
of the virial radius. In turn, we confirm that the baryon 
fraction apparently increases significantly from the center 
to outer parts of clusters. 

Several previous studies have shown a clear trend of in- 
creasing (3 with Tx, while we do not find such a clear 

trend, neither we find a trend of fh with /3 which is con- 
sistent with the absence of an fb-Tx correlation, although 
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we find a slightly decreasing fl^^ (derived from the IHE 
method) with increasing /3, but with a large dispersion. 
This result is important as it shows that the scatter in the 
baryon fractions derived from the IHE method is probably 
not due to the sole errors in the j3 measurement, but is 
rather due to the IHE mass estimator itself. 
Our sample does not show any evidence of the strong in- 
dication highlighted by D90 that in low temperature sys- 
tems, a larger fraction of baryons is present in the stellar 
component. Although a trend could be present in our sam- 
ple, the data are certainly consistent with a stellar to gas 
mass ratio being constant with temperature (or mass), a 
further argument that non gravitational processes such as 
galactic winds are playing a minor role in the overall dis- 
tribution of gas in clusters. 

Finally, it appears that the properties of X-ray clusters 
are still difficult to quantify because of the lack of large 
homogenous samples of clusters for which both optical and 
X-ray data are available. Such a situation is likely to im- 
prove with Chandra and XMM. Nevertheless, the sample 
we have studied reveals that clusters show important dif- 
ferences in the detail of the structure of their baryonic 
content, but that their global properties, baryonic frac- 
tion and stellar content do not show strong systematic 
differences with temperature. 

Notes on individual clusters : 

- A85 : Optical data for this cluster are unsafe and ex- 
tend only out to 900 kpc. We used the observations of 
Murphy (1984) but, as the fit with an unusual galaxy 
density profile he performs is rather poor (and suffers 
from an inconsistency between Hq = 50 and Hq = 60), 
we chose to replace it with a standard King profile. 

- A401 : As Buote & Canizares (1996) do not give the 
central electron density, we computed it with our pro- 
gram in the samc^ way as for Abell 665, since both 
objects were observed with ROSAT PSPC, using the 
galactic hydrogen column density of David et al. (1993) 
and the count rate inside a given radius provided by 
Ebeling et al. (1996). 

- A2029 : The same as for A401 applies. 

- A2163 : Optical data for this cluster are unsafe. No 
galaxy distribution was available. We therefore fitted 
the integrated luminosity profile given in Squires et al. 
(1997), but it was not corrected for background galax- 
ies and it extends only to 1.3 Mpc whereas i?xiim = 
4.6 Mpc. 

- AWM7 : We used the list of galactic positions and 
magnitudes within 1° of the central cD of Beers et al. 
(1984) to build an integrated luminosity profile, cor- 
rected for incompleteness using their limiting magni- 
tude and the standard Schechter luminosity function. 
The optical core radius was imposed to be the same as 
the X-ray core radius, which gives very similar results 
as excluding the three innermost galaxies (otherwise. 



the fitted core radius is too small and in fact, the King 
form is not a good representation of the central parts 
of clusters). 

- Hydra A : The same procedure as for A2163 was ap- 
plied, with the three points given by D90 : Ly / Lyq = 
8.2 10" at 0.5 Mpc, 1.3 lO^^ l Mpc and 1.9 lO^^ 

2 Mpc. 

- HCG 62 : We derived an integrated number count pro- 
file from the list of galactic positions of Zabludoff & 
A/fulchaey (1998) (using their velocity criteria to select 
true members) and fitted it with the function in Equ.7, 
assigning to each galaxy the mean luminosity derived 
from the limiting magnitude of the observations and 
the standard Schechter luminosity function. The num- 
ber of galaxies contained in this "compact group" is 
much larger than usually assumed (45 members with 
tob < 17 instead of 4 in Hickson 1982). 

- HCG 94 : EbeUng et al. (1995) claim this object has 
been misclassified and, from its X-ray emission, looks 
more like a poor cluster rather than a compact group. 
Only 7 galaxies are generally attributed to HCG 94 
but we made use of the indication of Ebeling et al. 
that 12 more galaxies are observed within a 1 Mpc 
radius and at me < 18, to which we attribute a mean 
luminosity as for HCG 62. The fit by Equ.7 we perform 
relies entirely on this point (fixing the core radius at 
the X-ray value) since inclusion of the central galaxies 
would lead to a physically unacceptable core radius). 
Therefore, optical data for this object are unsafe. 

- NGC 533 : As no central electron density is given by 
Mulchaey et al. (1996), we computed it from the gas 
mass that they obtain at a given distance. For the opti- 
cal part, this is the same case as HCG 62. This group 
contains 36 members with tob < 17 instead of 4 in 
Geller & Huchra (1983). 

- NGC 2300 : The same as for AWM7 applies, using 
magnitudes from the RC3 and excluding the two cen- 
tral galaxies from the fit instead of fixing the core ra- 
dius. 

- NGC 4261 : Since the central electron density Davis et 
al. (1995) give is inconsistent with their total gas mass, 
we computed it (this is again the same case as for A401 
and A2029). We also used optical data directly from 
Nolthenius (1993) and applied the same method as for 
AWM7 with magnitudes taken from the RC3 (except 
we did not have to impose the core radius). 

- RXJ 1340.6 has not been included in figures showing 
mass ratios as a function of overdensity, because it is 
a very peculiar case : the interior of the central giant 
elliptical galaxy is seen through a very large range of 
overdensities (at least out to S = 7000). 

- Several clusters have unreliable X-ray temperatures : 
A76, A426 (very strong coohng flow), A1377, A1775 
(likely very strong cooling flow) and A2218 (steeply 
outwards- decreasing temperature profile). 
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- Error bars on gas mass for groups include an estimate 
of the metallicity uncertainty, which results in an un- 
certainty on the electron density. This effect was taken 
into account only for groups, because it is significant 
mostly in the case of low temperatures. 
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Table 1. X-ray data. Parameters refer to Eq. 9 : rico is the central electron density. Temperature error bars of clusters 
are computed from different reliable references (taking into account or not a cooling flow) as the maximum of two 
estimates : the dispersion among the measures and the quadratic mean of the quoted uncertainties (this procedure 
produces large uncertainty when there exists a possibility for strong temperature gradient). They are given at a 90% 
confidence level, multiplying when necessary la errors by 1.64 (the errors of Fornax and RXJ1340.6 given by 196 and 
PA94, which confidence levels are not stated, are assumed to be la). 

*^*' Fornax and HCG62 have got a two-component gas density profile, the one in the second line corresponding to the 
nucleus. When no limiting radius i?xiim is given, we assumed i?xiim (Mpc) ~ Tx/3 (keV), a relation calibrated on 
other clusters. 



name 


Z 


Jx 

yKG V ) 


ref 




P 


fcX 

^^ivipc ) 


TleO 

^^lu cm J 


ref 


A76 


0.0416 


1 p:+2.8 


D93 




0.58 


0.356 


1.07 


CNT97 


A85 


0.0518 


6.1 ±0.44 


MF98 


1.4 


0.497 


0.06 


6.5 


P97 


A119 


0.0440 


5.62 ±0.61 


WOO 




0.56 


0.378 


1.18 


CNT97 


A401 


0.0748 


8.3 ± 1.16 


MF98 


2. 


0.606 


0.285 


4.31 


BC96, EVB96 


A426 (Perseus) 


0.0183 


6.79 ± 0.76 


F98 


2.5 


0.727 


0.416 


2.89 


CK97 


A576 


0.0381 


4.02 ± 0.31 


WOO 


1.5 


0.64 


0.24 


2.47 


MG96 


A665 


0.1816 


8.26 ± 1.29 


D93 


2.36 


0.763 


0.44 


2.85 


this study 


A1060 (Hydra I) 


0.0114 


3.1 ±0.14 


T96 




0.61 


0.094 


5.00 


LM96 


A 1377 


0.0514 


2.6«:^ 


W.IF97 




0.46 


0.188 


0.75 


CNT97 


A1413 


0.1427 


7.32 ± 1.10 


WOO 




0.62 


0.156 


8.00 


CNT97 


A1656 (Coma) 


0.0232 


8.11 ±0.79 


D93 


4. 


0.75 


0.42 


2.89 


BHB92 


A1689 


0.1810 


9.23 ± 1.05 


WOO 




0.72 


0.205 


14.0 


CNT97 


A 1775 


0.0696 


3.69 ± 3.03 


WOO 




0.58 


0.174 


2.73 


CNT97 


A2029 


0.0767 


8.7 ±0.75 


MF98 


2.3 


0.682 


0.251 


6.5 


BC96, EVB96 


A2052 


0.0348 


3.03 ± 0.23 


WOO 




0.67 


0.142 


4.43 


CNT97 


A2063 


0.0337 


3.52 ± 0.64 


woo 




0.50 


0.074 


5.9 


CNT97 


A2163 


0.201 


14.6 ± 1.21 


EAB95 


4.6 


0.62 


0.305 


6.65 


EAB95 


A2199 


0.0303 


4.27 ±0.28 


WOO 




0.62 


0.117 


6.5 


CNT97 


A2218 


0.175 


6.84 ± 2.83 


WOO 


2.1 


0.71 


0.276 


5.54 


SK96 


A2256 


0.058 


6.98 ± 0.44 


WOO 


3. 


0.795 


0.54 


2.51 


HBN93 


A2634 


0.0312 


3.27 ±0.28 


woo 




0.58 


0.32 


0.89 


CNT97 


A2657 


0.0414 


3.81 ±0.40 


woo 




0.52 


0.124 


4.4 


CNT97 


A2670 


0.0745 


3.73 ± 0.70 


woo 




0.70 


0.174 


3.83 


CNT97 


AWM7 


0.0176 


3.79 ± 0.36 


woo 


1.25 


0.53 


0.102 


6.47 


NB95 


Hydra A (A780) 


0.0522 


3.56 ± 0.37 


woo 


1. 


0.7 


0.145 


6.5 


D90 


Fornax'*-* 


0.0046 


1.09 ±0.18 


196 


0.4 


0.60 

0.51 


0.175 

0.007 


0.70 

19.6 


196 


HCG62(*) 


0.0138 


0.96 ± 0.07 


PB93 


0.36 


0.38 
0.9 


0.06 
0.019 


1.53 
15.6 


PB93 


HCG94 


0.04218 


2.75tr4 


EMW95 


1. 


0.58 


0.17 


3.0 


EMW95 


NGC533 


0.0172 


1 05+"-"" 


MD96 


0.58 


0.69 


0.237 


0.7 


MD96 


NGC2300 


0.0076 


0.97t'>:iJ 


DM96 


0.33 


0.41 


0.057 


3.11 


DM96 


NGC4261 


0.0088 


0.851°;?^ 


DM95 


0.6 


0.31 


0.026 


2.9 


DM95 


NGC5044 


0.0087 


0.98 ± 0.04 


D94 


0.4 


0.53 


0.028 


7.9 


D94 


RXJ 1340.6-F4018 


0.171 


0.92 ± 0.08 


PA94 


0.4 


1. 


0.181 


2.5 


PA94 
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Table 2. Optical data. The first part corresponds to the galaxy spatial distribution and the second part to the 
luminosity function (Eq. 5). 

(*) TOP : type of optical profile : (1) stands for a King form (Eq. 2), (2) for a de Vaucouleurs form (Eq. 3), (3) for an 

integrated luminosity profile (Eq. 7) and (4) for a dc Vaucouleurs projected luminosity density profile (cf Eq. 3). 
(*) ao : (1) and (2) :' galaxies Mpc^^ ; (3) : Lq (IO^Lb©) ; (4) : ctlo (10"Lbo Mpc""). 

Tc : for (2) and (4), corresponds to Vy. 

e : for (2) and (4), corresponds to 7. 



ri (\ irio 




''d 


,.(•) 

' c 

(Mpc) 




M-n 1- 
-''-^B lim 


rcf 






rcf 


A76 


(2) 


39.8 


0.485 


0.678 


-18.77 


CNT97 


-21.00 


1.25 


CNT97 


A85 


(1) 


253 


0.518 


1 


-17.76 


M84 


-21.17 


1.13 


OH89 


A119 


(2) 


214 


0.023 


0.271 


-19.77 


CNT97 


-20.98 


1.25 


CNT97 


A401 


(1) 


165 


0.4 


1 


-18.27 


D78b 


-20.96 


1.25 


D78a 


A426 


(1) 


209 


0.308 


1 


-17.34 


KS83 


-20.93 


1.25 




A576 


(1) 


250 


0.53 


1 


-16.26 


MG96 


-20.58 


1.18 


MG96 


A665 


(1) 


168 


0.5 


1 


-18.27 


D78b 


-20.21 


1.25 


D78a 


A1060 


(1) 


228 


0.180 


1 


-16.62 


B078 


-20.93 


1.25 




A1377 


(2) 


233 


0.170 


0.751 


-19.47 


CNT97 


-21.00 


1.25 


CNT97 


A1413 


(2) 


281 


0.05 


0.5 


-19.97 


CNT97 


-20.59 


1.25 


CNT97 


A1656 


(1) 


200 


0.34 


1 


-18.25 


KG82 


-20.93 


1.25 




A1689 


(2) 


140 


0.059 


0.43 


-20.77 


CNT97 


-21.66 


1.25 


CNT97 


A1775 


(2) 


205 


0.007 


0.227 


-19.57 


CNT97 


-21.44 


1.25 


CNT97 


A2029 


(1) 


198 


0.35 


1 


-18.27 


D78b 


-21.14 


1.17 


OH89 


A2052 


(2) 


175 


0.105 


0.553 


-19.57 


CNT97 


-20.35 


1.25 


CNT97 


A2063 


(2) 


153 


0.19 


0.767 


-19.47 


CNT97 


-20.30 


1.25 


CNT97 


A 01 do 


(3) 


1 on 


n coo 
U.OOZ 














A2199 


(2) 


442 


0.008 


0.264 


-19.27 


CNT97 


-20.17 


1.25 


CNT97 


A2218 


(1) 


387 


0.4 


1 


-18.27 


D78b 


-20.52 


1.25 


D78a 


A2256 


(1) 


138 


0.49 


1 


-18.27 


D78b 


-21.27 


1.12 


OHJ87 


A2634 


(2) 


508 


0.008 


0.252 


-19.27 


CNT97 


-19.75 


1.25 


CNT97 


A2657 


(2) 


648 


0.01 


0.345 


-19.57 


CNT97 


-20.86 


1.25 


CNT97 


A2670 


(2) 


245 


0.046 


0.393 


-19.47 


CNT97 


-21.00 


1.25 


CNT97 


AWM7 


(3) 


6.25 


0.102 






B84 








Hydra A 


(3) 


6.07 


0.146 






D90 








Fornax 


(1) 


250 


0.337 


1 


-13.20 


F89 


-20.82 


1.32 


FS88 


HCG62 


(3) 


3.07 


0.180 






ZM98 








HCG94 


(3) 


7.72 


0.17 






HKA89, EMW95 








NGC533 


(3) 


2.64 


0.140 






ZM98 








NGC2300 


(3) 


3.26 


0.245 






G93 








NGC4261 


(3) 


5.82 


0.163 






N93 








NGC5044 


(1) 


331 


0.188 


1 


-15.34 


FS90 


-20.93 


1.25 




RXJ 


(4) 


3449 


0.00524 


0.665 




PA94 









H. Roussel, R. Sadat, A. Blanchard: Baryon Content of Groups and Clusters 



19 



Table 3. Dynamical quantities for the whole sample, at the limiting radius r2oo and with the rsoo -Tx normalization 
of EMN (the redshift is taken into account in this relation). The scalings with the Hubble constant are : M^ot and 

-1 -- -2 

M* oc /150 , Mgas oc h^Q and L oc h^^ . 



name 


''200 


i\/r 

Mtot 


i\/r 




f 

Jh 


T\/f / T 
M 1 Ij-b 


jWihe/jWslm 




(MpC ) 


/"in 
(lU 


Mq) 




(/o) 






A76 


1.37 


20.12 


5.42 


0.27 


28.3 


239 


0.61 


A85 


2.72 


129.70 


10.21 


2.90 


10.1 


143 


0.70 


A119 


2.64 


125.44 


19.21 


1.97 


16.9 


204 


0.72 


A401 


3.07 


225.02 


44.32 


1.56 


20.4 


462 


0.75 


A426 


3.01 


177.79 


33.23 


0.81 


19.2 


700 


0.91 


A576 


2.25 


75.29 


10.74 


0.99 


15.6 


243 


0.84 


A665 


2.65 


183.81 


27.92 


1.60 


16.1 


367 


0.98 


A1060 


2.05 


49.42 


4.39 


0.24 


9.4 


663 


0.87 


A1377 


1.77 


37.30 


3.87 


0.96 


13.0 


124 


0.62 


A1413 


2.63 


158.22 


21.86 


1.01 


14.5 


499 


0.83 


A1656 


3.26 


221.40 


33.66 


1.46 


15.9 


486 


0.98 


A1689 


2.81 


223.32 


38.49 


2.43 


18.3 


294 


0.92 


A1775 


2.06 


62.65 


8.33 


0.83 


14.6 


242 


0.78 


A2029 


3.13 


230.00 


36.53 


1.72 


16.6 


427 


0.89 


A2052 


1.96 


48.05 


5.38 


0.84 


13.0 


183 


0.91 


A2063 


2.12 


61.68 


8.43 


0.79 


15.0 


250 


0.67 


A2163 


3.44 


443.64 


82.83 


6.97 


20.2 


204 


0.77 


A2199 


2.34 


80.81 


9.24 


1.13 


12.8 


228 


0.85 


A2218 


2.44 


144.83 


24.80 


2.78 


19.0 


167 


0.89 


A2256 


2.88 


179.35 


35.88 


2.02 


21.1 


285 


0.95 


A2634 


2.05 


55.71 


7.14 


1.43 


15.4 


125 


0.75 


A2657 


2.18 


71.43 


12.26 


0.76 


18.2 


303 


0.67 


A2670 


2.06 


61.11 


6.16 


0.96 


11.6 


204 


0.97 


AWM7 


2.25 


73.61 


13.05 


0.56 


18.5 


423 


0.68 


Hydra A 


2.07 


59.78 


7.38 


0.46 


13.1 


420 


0.95 


Forna:x 


1.23 


10.83 


1.29 


0.10 


12.8 


337 


0.78 


HCG62 


1.14 


9.34 


1.50 


0.15 


17.6 


196 


0.50 


HCG94 


1.85 


43.90 


7.61 


0.51 


18.5 


273 


0.74 


NGC533 


1.18 


10.01 


1.10 


0.15 


12.5 


207 


0.91 


NCG2300 


1.16 


10.33 


2.31 


0.13 


23.7 


249 


0.49 


NGC4261 


1.08 


8.68 


1.92 


0.30 


25.6 


93 


0.36 


NGC5044 


1.16 


8.98 


0.81 


0.17 


10.9 


171 


0.74 


RXJ 


0.90 


6.65 


0.69 


0.15 


12.6 


144 


1.32 
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Table 4. Average dynamical quantities for the objects with the most reliable data (details on those that have been 
discarded can be found in the paragraph "notes on individual clusters"), using the EMN normalization. Values are 
given virial radius rgoo (from SLM) but also at two other ones : r^oo, the limiting radius within which EMN claim 
that the hydrostatic equilibrium is universally reached ; r2ooo which we preferred to use because this represents the 
maximal extent of X-ray observations that is valid for the whole sample, including groups. Groups are defined by 
Tx < 2 keV and hot clusters by Tx > 5 keV. 
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DM profile) 
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at r2oo 


14.6 


16.0 


13.6 


93.3 
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20.1 


22.0 


77.3 


247 


all 


at rsoo 


11.9 


13.4 


10.5 


86.9 


278 


17.0 


19.2 


66.6 


213 




at r2ooo 


8.9 


10.4 


7.7 


84.5 
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12.8 


15.2 


60.9 


195 




at r2oo 


14.7 


16.5 


9.0 


62.4 


200 


26.6 


29.8 


45.7 


146 


groups 


at rsoo 


11.0 


13.1 


6.3 


58.3 


186 


19.7 


23.4 


39.0 


125 




at r2ooo 


7.0 


9.6 


4.2 


61.1 


196 


12.3 


16.6 


37.5 


120 




at r2oo 


13.4 


14.6 


12.2 


95.0 


304 


17.2 


18.8 


77.8 


249 


cool clusters 


at rsoo 


11.0 


12.3 


9.3 


88.6 


283 


15.0 


16.7 


67.9 


217 




at r2ooo 


8.4 


9.7 


6.9 


84.6 


271 


11.6 


13.5 


63.2 


202 




at r2oo 


15.9 


17.0 


19.4 


118.1 


378 


18.8 


20.2 


104.3 


334 


hot clusters 


at rsoo 


13.7 


14.8 


15.7 


109.8 


351 


17.3 


18.9 


89.1 


285 




at r2ooo 


10.7 


11.9 


11.8 


104.8 


335 


14.4 


16.0 


78.5 
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Table 5. Same as Table 3 but with the r2oo-^x normalization of BN. 
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name 


T'200 


Mtot 


Mgas 




/b 


M/Lb 


Mihe/Mslm 




(Mpc) 


(10 


"Mo) 




(%) 


(Mq/Lbq) 




A 76 


1.52 


26.47 


6.40 


0.31 


25.4 


272 


0.52 


A85 


3.02 


174.74 


11.97 


3.11 


8.6 


180 


0.58 


A119 


2.93 


167.47 


22.42 


2.27 


14.7 


236 


0.60 


A401 


3.41 


297.72 


50.75 


1.65 


17.6 


577 


0.63 


A426 


3.34 


234.89 


37.20 


0.86 


16.2 


879 


0.77 


A576 


2.50 


100.32 


12.21 


1.07 


13.2 


299 


0.70 


A665 


2.95 


244.16 


31.16 


1.72 


13.5 


455 


0.83 


A1060 


2.28 


66.58 


4.99 


0.25 


7.9 


852 


0.72 


A1377 


1.97 


50.07 


4.63 


0.99 


11.2 


162 


0.52 


A1413 


2.92 


211.20 


24.82 


1.04 


12.2 


650 


0.69 


A1656 


3.62 


294.04 


37.42 


1.53 


13.2 


615 


0.82 


A 1689 


3.12 


295.02 


42.65 


2.61 


15.3 


362 


0.78 


A1775 


2.29 


83.81 


9.59 


0.96 


12.6 


278 


0.64 


A2029 


3.48 


305.31 


40.92 


1.82 


14.0 


537 


0.75 


A2052 


2.18 


64.21 


6.04 


0.89 


10.8 


230 


0.76 


A2063 


2.35 


82.53 


9.89 


0.81 


13.0 


328 


0.56 


A2163 


3.82 


586 39 


94.45 


7.42 


17.4 


253 


0.65 


A2199 


2.60 


108.21 


10.48 


1.28 


10.9 


270 


0.70 


A2218 


2.71 


191.28 


27.75 


2.97 


16.1 


206 


0.75 


A2256 


3.20 


235.80 


39.93 


2.15 


17.8 


350 


0.81 


A2634 


2.28 


74.50 


8.31 


1.64 


13.4 


146 


0.63 


A2657 


2.42 


95.12 


14.32 


0.80 


15.9 


378 


0.56 


A2670 


2.28 


81.87 


6.87 


1.06 


9.7 


248 


0.80 


AWM7 


2.50 


97.92 


15.18 


0.58 


16.1 


542 


0.57 


Hydra A 


2.30 


79.83 


8.21 


0.48 


10.9 


537 


0.79 


Fornax 


1.37 


14.52 


1.49 


0.11 


11.0 


414 


0.65 


HCG62 


1.26 


12.51 


1.82 


0.16 


15.8 


246 


0.41 


HCG94 


2.05 


58.31 


8.77 


0.54 


16.0 


345 


0.62 


NGC533 


1.32 


13.41 


1.25 


0.16 


10.6 


263 


0.76 


NCG2300 


1.28 


13.73 


2.79 


0.14 


21.4 


308 


0.41 


NGC4261 


1.20 


11.55 


2.39 


0.32 


23.4 


117 


0.30 


NGC5044 


1.29 


12.07 


0.94 


0.18 


9.3 


216 


0.61 


RXJ 


1.00 


8.85 


0.74 


0.15 


10.0 


194 


1.11 
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Table 6. Same as Table 4, but using the BN normalization. 
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28.8 


31.8 


48.2 


154 


groups 


at rsoo 


9.7 


11.4 


6.8 


72.1 


231 


21.2 


24.8 


40.3 


129 




at r2ooo 


6.3 


8.3 


4.4 


71.7 


230 


13.2 


17.4 


36.8 


118 




at r2()o 


11.5 


12.5 


13.1 


119.7 


383 


17.7 


19.3 


81.4 


260 


cool clusters 


at rsoo 


9.4 


10.5 


9.9 


110.0 


352 


15.5 


17.3 


69.9 


224 




at r2ooo 


7.3 


8.4 


7.1 


101.7 


325 


12.1 


14.0 


63.1 


202 




at r2oo 


13.5 


14.4 


20.6 


147.7 


473 


19.1 


20.5 


109.5 


350 


hot clusters 


at rsoo 


11.7 


12.6 


16.5 


135.5 


434 


17.7 


19.2 


92.2 


295 
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9.3 


10.3 


12.3 


125.2 


401 


14.9 


16.5 


79.2 


253 



